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ABSTRACT 

There are many flows of practical importance where both Tollmien- 
Schlichting waves and Taylor-Gortler vortices are possible causes of tran- 
sition to turbulence. In this paper, the effect of fully nonlinear Taylor- 
Gortler vortices on the growth of small amplitude Tollmien-Schlichting waves 
is investigated. The basic state considered is the fully developed flow 
between concentric cylinders driven by an azimuthal pressure gradient. It is 
hoped that an investigation of this problem will shed light on the more com- 
plicated external boundary layer problem where again both modes of instability 
exist in the presence of concave curvature. The type of Tollmien-Schlichting 
waves considered have the asymptotic structure of lower branch modes of plane 
Poisseulle flow. Whilst instabilities at lower Reynolds number are possible, 
the latter modes are simpler to analyze and more relevant to the boundary 
layer problem. The effect of fully nonlinear Taylor-Gortler vortices on both 
two-dimensional and three-dimensional waves is determined. It is shown that, 
whilst the maximum growth as a function of frequency is not greatly affected, 
there is a large destabilizing effect over a large range of frequencies. 
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1. INTRODUCTION 


In laminar boundary-layer flows over a surface, such as a wing, shear- 
flow InstablllCles in the form of Tollmien-Schlichting waves can occur. These 
waves are the subject of much theoretical and experimental Interest since it 
is thought that they cause transition to turbulence. When the flow is over a 
curved surface, centrifugal instabilities such as Taylor or Gortler vortices 
may also be present. These may interfere destructively with the Tollmien- 
Schlichting waves and thereby delay transition. Alternatively, by making the 
flow three-dimensional, they could play an essential part in the process of 
transition. The interaction of these two types of instabilities is therefore 
of some theoretical importance and has practical applications in the develop- 
ment of laminar-flow wings. 

Hall and Bennett (1986) showed that when Tollmien-Schlichting waves 

travel past a curved boundary, an unstable Stokes-layer forms on the wall, and 

it was suggested there that the growth of longitudinal vortices in this 

Stokes-layer could destroy the Tollmien-Schlichting waves. In this paper, we 

consider the opposite problem, namely the stability of a Dean (1928) type 

Taylor vortex in a channel to small amplitude travelling waves. By comparing 

our results with the stability analysis of a channel flow without any vortex 

motion, we hope to be able to tell whether the presence of the vortices 

hinders or enhances the growth of the Tollmien-Schlichting waves. A related 

• • 

problem was studied by Nayfeh (1981) in which Gortler vortices were allowed 
to Interact with oblique Tollmien-Schlichting waves. There the Gortler vortex 
was determined by solving the parallel flow linear instability equations and 
had its amplitude assigned arbitrarily. Such a procedure could lead to 
Incorrect results because Hall (1982a, 1983) has shown that nonparallel 
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effects cannot be ignored in the linear Gortler instability problem. 
Moreover, a finite amplitude Gortler vortex has its amplitude determined by 
the Gortler number and cannot be specified arbitrarily. Furthermore, in the 
only case where a nonparallel theory of nonlinear Gortler vortices has been 
given (Hall (1982b)), the mean flow distortion induced by the fundamental is 
the same size as the fundamental. In such a situation, it is clear that the 
contribution of the mean flow correction and its harmonics cannot be 
ignored. The channel flow considered here does not vary in the streamwise 
direction, so non-parallel effects do not occur. Bennett (1986), however, has 
shown that our analysis does apply to external non-parallel flows, though 
results for that problem will not be available until the fully nonlinear 
Gortler problem in external flows has been solved numerically. 

We confine our attention to the linear stability of the vortex motion at 
high Reynolds numbers. Furthermore, we shall concentrate on the lower branch 
of the neutral curve, so that the Tollmien-Schlichting waves are governed by 
interactive boundary-layer theory. This case describes asymptotically almost 
all the unstable range of high Reynolds number disturbances. In Section 2, we 
derive the dispersion relation linking the wave-frequency to the wave-number 
for waves travelling parallel to the main direction of the flow. This is done 
in a similar manner to Smith (1979a), where the stability of unidirectional 
flow was considered. The difference between that and the present work is that 
there the basic flow varied on a cross-stream (z) distance comparable to the 
long wavelength of the disturbances, whereas the z-variation in our basic flow 
is governed by the shape of the Taylor vortex. In Sections 2 - 6, we consider 
"square" Taylor vortices, where the z variation is comparable to the channel 
width and therefore much faster than the streamwise (x) variation of the 
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waves. In Sections 3 and 4, we look at two limits of the dispersion relation 
derived in Section 2. Firstly, in Section 3, we look at what happens when the 
amplitude of the Taylor vortex is small, so that the vortex is governed by the 
weakly non-linear theory of Seminara (1976). Secondly, in Section 4, we find 
how the scaled wave number a of the Tollmien-Schlichting waves behaves 
when the scaled wave frequency is large. In Section 5, we describe the 
numerical calculations needed to work out the vortex velocity field and to 
find solutions of the dispersion relation. Section 6 extends these results to 
the case of waves travelling obliquely to the flow. Finally, in Section 7, we 
give a discussion of our results and their relevance. 


2. THE DISPERSION RELATION FOR SHALL AMPLITUDE TOLLMIEN-SCHLICHTING HAVES IN 
THE PRESENCE OF FULLY NONLINEAR TAYLOR-GORTLER VORTICES 

We take as our basic flow the Taylor vortex that arises in the Dean 
(1928) problem when incompressible fluid is driven between concentric 
cylinders by a constant azimuthal pressure gradient. If the radii of the 
cylinders are a and a + d, then we assume that the channel is narrow, that 
is 6 = a/d » 1, so that the Taylor vortex is an instability of plane 

A ^ ^ 

Poiseuille flow, driven by centrifugal forces. If (r , 0 , z ) are 
cylindrical polar coordinates, with r* = 0 corresponding to the axes of the 
cylinders, we define dimensionless coordinates (x,y,z,t) by 

“ft ^ ^ ^ 

X = a0 /d, y = (r - a)/d, z = z /d, t = U^t /dRe 


where the Reynolds number Re = Ujj^d/v , v is the viscosity, is a typical 
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raean flow speed, and t is a dimensionless time. The dimensionless velocity 
and pressure of the Taylor vortex, (u,v,w) and p, are given by 


u* = U (U„ + u, 

— ra 0 2Re 


V, ^w). P 


* .2. 12x . p V 


( 2 . 1 ) 


and 


Uq = 6y(l - y) 


( 2 . 2 ) 


is the mean flow driven by the pressure gradient. 

Substituting these expressions into the Navier-Stokes equations and 

2 

ignoring terms of 0(6 ) and 0(Re6 ) , we get 


3 V ^ 

3y 3 z 


= 0 , 




0 


Nu, 


(V 


2 


3 no2 


+ TU 


0 




Nv + 



3y3z 


.isiu 


(2.3) 


where 



N 


- 1 /■ 3 u. 

= 2 S7 




(2.4) 


The boundary conditions are 


u = v = w = 0 on y = 0,l 


(2.5) 


and u, V, and w periodic in z. 


whilst the Taylor number T has been defined by 
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T = 4Re6^. (2.6) 

The linear Instability problem discussed by Dean (1928) can be obtained 
by linearizing (2.3), so that the right hand sides vanish, and by replacing 
9/9t by 0. If the Taylor number T is plotted against the wave-number 
k for steady solutions of Dean's problem, periodic in z with period 2iT/k, 
an open neutral curve typical of convective or centrifugal instabilities is 
found. Points in (k,T) space above this neutral curve correspond to 
unstable linear Taylor vortices, whilst those below represent Taylor vortices 
that decay to zero when t u . The critical point of the curve is given 

by T = Tj, = 5161.86, k = k^, = 3.951. Here we are Interested in fully 

nonlinear steady solutions of (2.2) - (2.6). These exist in a region above 
the linear neutral curve and are obtained numerically in the manner described 
in Section 5. 

We now consider what happens when the Taylor vortex velocity (2.1) is 

perturbed by high Reynolds number Tollmien-Schlichting waves travelling 
parallel to the x-axis. For Re » 1, the components of velocity in (2.1) 
perpendicular to the x-axis become negligible. Smith analyzed the stability 
of a unidirectional flow depending on two spatial variables y and z near 

the lower branch of the neutral curve when the perturbations vary on a slow 
X length scale of 0(Re^^^). In his work, the variation of the basic flow in 
the cross stream direction z was also on a long length scale of 0(Re^^^). 
In this section, the basic flow varies on a relatively fast 0(1) lengthscale 
in z, forced by the behavior of the Taylor vortex. There are, however, cir- 
cumstances in which the z variation of the Taylor vortex is of 0(Re^^^), 
but these occur at much higher Taylor numbers and therefore are not discussed 


here 
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Following Smith, then, but taking into account the different z scales, 
we write 

e = X = e~^X, T = U t*/d. (2.7) 

m 


The flow splits up into three regions, an Inviscid core, and a viscous 

2 

critical layer of thickness 0(e ) on each wall. 

In the core, we perturb the Dean problem as follows 


^ = Um[U,o,o) + (e^ u, v, w)E + ...] 

* 2 4“ 

p = pu;-(e p E + ...) 
m 


( 2 . 8 ) 


where U(y,z) = Uq + u is the velocity of the mean flow and Taylor vortex. 


E - h exp(i(aX - ^^T)), h « 1 


(2.9) 


and the variables denoted by ^ are functions of y and z only. On sub- 
stituting into the Navier-Stokes equations, we get 


iau + v + w = 0 
y z 


iaUu + vU + wU = 0 

y 2 


i a U V = - p 

y 


i a U w = - p 


(2.10a-d) 


I 


1 
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with slipping conditions at the walls 


V = 0 at y = 0,1 


( 2 . 11 ) 


A A 




and u, V, w, and p are periodic in z. The scalings for u, v, and 

A A 

p are those of Smith whilst that for w is forced by a comparison of the 

last two momentum equations (2.10c,d). The velocities in (2.9) can now be 

written in terms of the pressure 


U. 


U 


(ia) u«[7 p-p— -p — ]/U 


y u ^ u 


A 


iov = - p /U 

y 


( 2 . 12 ) 


law = - Pjj/U 


Substituting these expressions into the first momentum equation (2.10b), we 
get 


o A P U p U 

,2 p. 2^- 2^-0. (2.13) 

u u 


The boundary conditions (2.11) together with the fact that U vanishes at 
both walls imply that 


Py = Pyy = = 1 ^t y = 0,1. (2.14) 


It can be shown from (2,13) and (2.14) that the core problem does not 
specify p uniquely, since for any solution of (2.13) and (2.14) we can add 
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on and multiply by arbitrary constants to get another solution. Thus p is 

determined by the interaction between the core and the viscous layers at y = 

0 , 1 . 

If p = <() is a solution of (2.13) and satisfies the boundary 
conditions 

4> = 0 at y = 0, (|) = I at y = 1, (2.15) 

A 

then it can be shown by series expansions in y and (1 - y) that p = <() 
also satisfies the boundary conditions (2.14). Hence 


P = Pq + (pQ ■ 

A 

is also a solution of (2.13) and (2.14) for arbitrary constants p^ and 
p^ and is therefore the general solution. 

As we move into the lower boundary-layer, the core pressure and unper- 
turbed velocity are such that 

p Pq + 0(y ), U -► XQ(z)y, as y + 0. (2.16) 

Hence from (2.12) the disturbance velocities are such that 


u 


Aq(z), 


-ioA^y, w^O(y^) as y 0. 


(2.17) 


The displacement term Aq satisfies 

X0(i«)^AQ = I (P^- Pq) <j)yyy|Q . 


(2.18) 
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The core velocities and pressure behave in a similar manner to (2.16) - (2.18) 
as we move into the upper boundary layer. 

In the lower critical layer, we write 


y = zh 


^ = U^[ (e ^XqY, 0,0) + (e^ u, v, w)E...] 


* 2 4 ~ 6 ~ 

p = pU^ [(e Po + e Pq^ + ...]E 


m 


O 

which, ignoring terms of O(h^), leads to the equations 


/V 


iau + v„ + w =0 
Y z 


rsf fst /St tst 


i(-n + uXqY)u + XqV + Aq^Yw = -iapQ + u^y 


tSJ fst 


i(-£2 + uXqY)w = -Pqj^ + w^Y 


(2.19a-e) 


Pq = const 


rK0 rsj 


Poi • 


with boundary conditions 


rst /N/ ^ 


u = v“W=0 on Y = 0, 


tsj ^ rst 


u + Aq, w -► 0, Pq + Pq as Y > «, 


(2.20a-c) 


rst fsj rst 


U. V, w, Pq, and p^^ periodic in z, 


- 10 - 


The scalings of u, v, and come from matching with the core 

whereas the scalings for w and have been chosen as large as 

possible. These two terms will be driven by matching with higher order terms 
In the core and are not specified uniquely by (2.19) and (2.20). Equations 
(2.18) differ from the linearized boundary-layer equations solved by Smith 
only in that the pressures p^ and p^^ are not equal here. Following 

Smith, (2.19c) can be solved for w in terms of Pqi* 

w = Pq^ KU)/aaXQr^ (2.21) 

z 


where 


M(5) = Ai(?) 


? s 

/ — ^ / Ai(s )ds 

Ai (s) « 


K = (iaXQ)^^^(y - Q/aX^), 


( 2 . 22 ) 


and Ai is the Airy function that satisfies Ai" = ^Ai. Eliminating v 
between (2.19a) and (2.18b) by differentiating (2.18b), substituting for w 
from (2.21) and solving for u we get 


(iaX^)^/^ u = -XqPqj M'(?) +^Poi {3M'(?) 

zz z 


2Ai'(0[~^|^]} + BAi(0 


(2.23) 


where the constant B is determined by the outer boundary condition on u 
in (2.20), 
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B = (iaXp) 


5/3 ^0 


X- 

_ Oz 0 

~ PqI ^0 K(5.) ’ 


(2.24) 


K(5 ) = / Ai(5)d5. 

5n 


From matching with the core and putting Y = 0 in (2.19b) and 

using (2.23) and (2.24), we get an expression for p^ in terras of p^^ 
and Aq. Substituting for Aq from (2.18) gives 


* Voi^ - ^ - '^0>Nyyly-0 


(2.25) 


’**0 Xq ^2 2Ai(?Q) ^ 


The problem in the upper critical layer at y ® 1 is the same as that 

at y = 0 with (X^ (z), A^, p^, p^^, p^) instead of (Xq(z), Aq, Pq, 

Pqi, Pq) in equations (2.19) to (2.24), and with (2.18) replaced by 


X^(ia)^Aj 


= -i(Pl - 


P0^'*’yyyly=l* 


(2.26) 


This leads to 


P, -I <Po - PlHyyyI. 

n - - xf 4 ^ ^ 


Ai'(5i> 


(2.27) 


,1/3^,. , .2/3 

5 — — i Q/(aXj^) • 
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In the limit as 8 /3 z > 0(Re Pq^> would tend to 

3 cid so that (2«25) and (2.27) would become the coupled second order 

differential equations of Smith. Our equations as they stand are easier to 
deal with as they are only first order differential equations in p^^ and 

~ z 

p^p Integrating (2.25) once and using the periodicity condition (2.20c) 
gives 


z 



A A 1 

+ (P, - Pq^ j/ 


2tt /k 


H<o(Zi)dZi 




yyy'O 


(IcOqT^^K^ 


dz , 


Integrating (2.27) gives a similar equation for p^ and p^ involving 

A A 

Cp and on eliminating p^ and p^ from thse two equations we obtain 


2n /k 


2a" = 


I 

0 




Ai'(?o) 


(iaXQ)^^\( 5 Q) 


yyy 


n dz/ 
y=0 


2n /k 


/ 

0 


V"l 


dz 


z 

2Tr/k ^ Ai'(5 ) 

0 (iax/73^^yyly=l 


z 



1 

dz. 


(2.28) 


This is the dispersion relation or eigenrelation, giving a in terms of 
fl for high Reynolds number, linear Tollmien-Schlichting waves. It applies 
to any unidirectional flow (U(y ,z) ,0,0) , so long as the period of the z 
variation 2ir/k « Re^^^. Thus the eigenrelation is also applicable to high 
Reynolds number flows in pipes of finite cross-section. Here z would 
correspond to distance measured around the pipe. 
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Three stages are needed to work out values of on the curve 

given by (2.28) for the Taylor vortex case U = Uq + u; 

(i) First, for a given Taylor number T we must find a steady solu- 
tion of (2.2) - (2.4) and hence U; 

(il) Once U is known we can determine the core pressure ((> by 
solving (2.13) and (2.15); 

(iii) Finally we can solve (2.28) with and given by 

(2.22), (2.25), and (2.27) and Xq and given by 


3. WEAKLY NONLINEAR THEORY 

In order to find nonlinear solutions of the vortex equations (2.2) - 

(2.5) , we need to use a numerical method, as in Section 5. However, when the 
Taylor number is only slightly greater than the critical linear Taylor 
number T^, the amplitude of the vortices is small and solutions of (2.2) - 

(2.5) are described by the weakly nonlinear theory of Seminara (1976). We now 
apply our dispersion relation (2.28) to Seminara"s velocity profile. The 
results we obtain indicate how the stability of the flow is affected when the 
flow becomes slightly three-dimensional and will provide a useful check on the 
full nonlinear calculations of Section 5. The weakly nonlinear velocity in 
the streamwise direction is given by 

— 2 3 

U = Uq + AUj^coskz + A (u2q + U22Coszkz) + 0(A ). 


(3.1) 
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Here the vortex wavenumber k is the critical one for linear vortices, k = 

= 3.951, and the vortex amplitude A is related to the Taylor number by 

T - T 1/2 

A = 0.1725( « 1. (3.2) 

c 

The velocities uj, ^22 given by Serainara (1976) and are 

independent of z . 

The solution of the pressure equation (2.13) with boundary conditions 
(2.15) is then forced to behave in a similar manner to (3.1) 

(|) = <|)q + A(J)^coskz + (3.3) 


Here we are only interested in the leading order effects of the vortex on the 

dispersion relation. Since the integral of coskz over a period is zero, the 

fundamental 0(A) terms will only appear as a product with other 0(A) 

2 

terms and so will only have an 0(A ) effect. Hence we must also take into 

account the mean flow correction terms U 2 q and ^20* ^ similar argu- 

2 3 

raent the first harmonic terms 0(A )cos2kz will only have an 0(A ) 

effect, and so are ignored here. Substituting (3.3) into (2.13) we find that 

the pressure term (J> is given by 


where 


5 4 3 

= 6y^ - 15y^ + lOy*" 


^0^20 ^1 '^l 5 2 

' 20 ^- 


1 U u - u, u, 

C./ 


(3.4) 


3 
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and that the fundamental pressure term is governed by 


(^)y 

“5 


- k 


2 


5 /l 


~2 “ 3 at y = 0,l, 


U, 


(3.5) 


If we write the skin friction at the walls as 


___ 2 
^0 “ ^yly-0 “ ^ ^UjCOskz + A (y2Q + •••) 

^1 “ "^yly^O “ ^ ^'’i^oskz + A^(v2q + ...), 


(3.6a, b) 


then from (3.4) the pressure terms we need to evaluate the dispersion relation 
are 


^0 !y=0 “ *0 ly=l “ 

yyy yyy 

^20 ly=0 ” ^^20 ~ 2^1 ‘•*1 ly=0 6^ ” 


(3.7) 


yyy 


yyy 


'*’20 ly=l “ ^^20 2 '’l '*’1 ly=l ^ 

yyy yyy ^ 


We can define a mean value for 5 q» 5^ by writing 




-Oe 


w/6 


(6a) 


273 


(3.8) 


and substituting (3.6a) into the definition of i|»q, (2.25), we find after 

some simplification that the integrating factor in (2.28) becomes 

z 


~ — (1 ~ ^ coskz + 0(A^)) 


(3.9) 
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where 

G = 3 + — ^ (?K + Al'). 
Ai(5) 


(3.9b) 


Also, the rest of the integrand in the eigenrelation can be written as 


Ai'CCg) 


1 


Au 


K( 5 Q)(ioXQ) 


T73 


77^ — F coskz + A^<7^ F, 

(6ia)^'^ K(C) ^ ^ 


(3.10) 


*^20 

18 


Fj)) 


where 


2 (5K + Ai") + 1 

K(5)A1'(5) 


(3.10b) 


and 


9F =1+5^ 

2 ^ Al' 


+ 1 «K. Ai-)(5 2« M 


-). 


(3.10c) 


Combining these results we can evaluate the part of (2.28) corresponding to 


the lower boundary layer 
z 

2Tr/k '^O'^^l Ai'(C^) —'‘'n 

/ e fTo <l> I 

i 1 /■? ryyy|y_Q J 


2n/k( V^l 


0 


(iox/^c^) 


0 


Ai'(0 60 n + A2r''l fp . 3. ^ 20/1 

^ (1 + A [_ (F^ + - 2> - — (— 


K(C) (6ia) 


- 2 ) 


h ly=0 

yyy 


Using a similar result for the other half of (2.28), we get the following 
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eigenrelation 


... vl/3 2 .—V 2 F-G « 

° - 1 - I- (F, f ‘ - |) ^ 

60Ai'(D ^ 2 12 2 

r 

(Mi4>1 Io+'’i'I>i li)(l - -^ - J)/360 - 2C) 


yyy 


yyy 


(3.11) 


where we have used the result ^20 ^ ^20 " constants In (3.11) are 

obtained from Seminara (1976) and by solving (3.5) 


li 


2 

1 



36 


8.072 


^^l^lyyylo ^ 1 ^ 1 yyy ^ 1 ^ 


6.843 


-2C = 26.86. 


We now choose to look only at the spatial stability of the vortex motion, 
that is for a given real frequency we solve (3.11) numerically to find 

the wavenumber. If a^, the imaginary part of the wavenumber, turns out to 
be positive, then the resulting waves decay downstream as x ->• », whilst 
if < 0 the waves will grow exponentially. Equations (3.8), (3.9b), 

(3.10b,c), and (3.11) were solved at various values of Q by Newton-Raphson 
iteration in a. Figure 1 shows plotted against for 

A = 0, .1, .15 corresponding to Taylor numbers of T = T^, 1.33T^, 1.76T^. 
The general pattern is the same in all three cases. For the frequency less 
than some critical value (depending on the Taylor number), all disturb- 
ances decay. At Q = ® that linear disturbances neither grow 
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nor decay. This polat corresponds to the asymptotic limit of the lower branch 
of the neutral curve. For all waves grow, and as U * <», (a^) 

decays to zero since we are tending towards the upper branch of neutral 
curve. Disturbances corresponding to the upper branch occur on different 
length and time scales to the lower branch disturbances, so that in an 
analysis near the upper branch x and t would be scaled on different powers 
of the Reynolds number to (2.7), so that however large is we will never 

actually reach the upper branch where = 0. The results of Figure 1 show 

that the vortices have negligible effect on the neutral frequency and 

that the growth rate over a large band of frequencies is significantly 
increased. This corresponds to a destabilization of Poisseulle flow by the 
vortices. We shall see that this trend is also found when the fully nonlinear 
problem is solved numerically. 


4. THE HIGH FREQUENCY LIMIT, Q -h • 

In this section we determine the asymptotic behavior of the eigenrelation 

(2.28) when ->■ «>. We assume that G/ ® as •♦• «, which means 

a 

that 5 q and ->■ «>. This can be checked at the end of the calculation. 
First we need two results giving the behavior of the Airy function, its deriv- 
ative, and its Integral for large arguments. Writing Ai(s) = Ai"/s and 
Integrating by parts, we find 


K(C) = / Al(s)ds 

5 


Al^(g) 

C 




(4.1) 


as C 


Also it is known (see e*g*, Abramowitz and Stegun (1965)) that the ratio of 
the Airy function to its derivative is given by 


Ai/Ai' as C •> 


(4.2) 


Using these two results we can calculate the asymptotic behavior of each term 
in (2.28). So from (2.25) using (4.1) and (4.2) we get 


Oz 1 


and using (2.22) we can write this as 


X^ a 

Oz Oz 




Hence the integrating factor 

z 

/ 'I' 0^2 

r\ 


0 


1 


XqU 


( 1 “ o 1 /o • • ♦ ) 5 


^0 (-ifi ) 


so that 


^dz 


2j>/k I *0“ 

/ e“ dk 


0 ‘o 


2n 

IT 




as 5 


Also from (4.1) and (4.2) 


Ai'(^0> 


X a 


(4.3) 


(4.4) 
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Corabining this with (4.3) gives 


e Ai (Cq) ^ ^ 

K(?Q)(iaXQ)^^^ oXq 


(4.6) 


Substituting (4.4) and (4.6) together with similar expressions involving C, 
into the dispersion relation (2.28), we obtain 


where 



2ira 


k(-i03)l/2 


L, - 


2Tra 




2iT/k 6 

-f ^1 

0 ^i 


dz 


y=i 


= / 
0 


2Tr/k 



(4.7) 


(4.8) 


We can now obtain an as}rmptotic series for a in inverse powers of 
from (4.7) in the form 




1/3 




-5/6 


(4.9) 


where 


$ $ 

<=? + '=2 


TT 

3kc^ 



(We note here that Q •> » implies that ^0* ^ assumed earlier.) 

Thus the imaginary part of a decays to zero like some constant 
times n The actual value of this constant depends on the vortex pro- 
file through the constants and C2* We now determine these constants for 
vortices governed by the weakly nonlinear theory described in Section 3. If 
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the velocity is given by (3.1), we find that 


and 


2tt 

6k 


(1 


2 ^1 

+ A - 



with a similar formula for L 2 but with v replacing p. 

2 2 

- .2 , y, + V, 

C3 = 10(1 (- i - ^ - i --) + 12c}) 


Hence 


^ 2 2 

20 

S ■ I7 » - r ^ -V-^» 


SO that the asymptotic form for the growth rate is given by 


~ /2 (10^^^)(l - l^C + I ( ^^35 


(4.10) 


This asymptote is plotted in Figure 1 for various values of A along 
with the corresponding weakly nonlinear dispersion relations. The asymptotic 
form (4.10) is seen to accurately predict over a wide range of fre- 
quencies and thus provides a useful check on the calculations of Section 3. 


5. THE NUMERICAL CALCULATION OF A FINITE AMPLITUDE TAYLOR VORTEX 

Here we describe how we integrated (2.3) numerically to find the finite 
amplitude Taylor vortex whose instability we wish to determine. The method 
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used is essentially that described by Rogers and Beard (1969) who investigated 
numerically the classical Taylor problem driven by the motion of the inner 
cylinder. Rogers and Beard solved a system similar to (2.3) by Fourier 
expanding u and v in the z direction and using finite differences in the 
radial direction. Later Fasel and Booz (1984) performed related calculations 
using finite differences in both directions. The method of the latter authors 
is apparently the most efficient at very high Taylor numbers where jet-like 
structures develop along the cylinders. Here we do not perform calculations 
at such high Taylor numbers, so we use the method of Rogers and Beard. 

Thus the velocity components in (2.3) are expanded as 


u = u^(y,t) cosknz, 


V = v„ + y V (y,k) cosknz, (5.1) 

U 1 n 


w = ][ w^(y,t) slnknz. 


Here we have anticipated the usual result that the only mean flow generated by 
the vortex is in the azumuthal direction. The expansions (5.1) are then 
substituted into (2.3), and the coefficients of cosknz are equated to give 
an infinite sequence of coupled nonlinear differential equations for 
{u^} and Thus for example the equation for Vq is 





_1 

2 


j=l 




(5.2) 


We obtain steady state solutions of (2.3) by integrating forward in time from 
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some appropriate initial guess. At Taylor numbers close to the critical, the 

initial guess can be taken to be the weakly nonlinear state described in 

Section 4. At high Taylor numbers the initial guess was taken to be the 

equilibrated solution from a previous calculation at a lower Taylor number. 

A fully implicit scheme was used to march forward in time. Hence if 

At is the time step and h the step length in the y direction we obtain a 

stable scheme for At ~ 0(h). The nonlinear terms on the right hand side of 

the disturbance equations were always evaluated explicitly. The number of 

Fourier modes and intervals in the y direction required to achieve a solu- 

tion sufficiently accurate for our purpose depends on the ratio ^ . In 

c 

order to monitor the energy in different harmonics, we followed Rogers and 
Beard and defined 


The number of axial modes required was varied until the converged values of 

Eq, Fq achieved sufficient accuracy. Similarly the step length h was 

varied until Eq, Fq converged to sufficiently accurate values. For the 

calculations reported here, it was found that eight axial modes and h = .01 

were sufficient to enable us to determine the dispersion relation to the 

accuracy indicated in Section 7. In Figures 2 and 3, we have shown the 

dependence of Eq and Fq on T for k = k^, = 3.951 the most dangerous 

Taylor vortex mode. We see that at sufficiently small values of - 1} 

c 

the results are consistent with the weakly nonlinear results which can be 
derived from Seminars (1976). However the fully nonlinear solution diverges 

m 

from the asymptotic result at quite small values of - l} so that the 

c 
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results predicted in Section 3 for A = .15 are beyond the range of 
validity of the Taylor vortex amplitude expansion. However we see in Figure 1 
that the full solution for T = 1.78 gives values of qualita- 
tively similar to the weakly nonlinear results for most values of the 
frequency shown in that figure. 

Once the Taylor vortex has been computed, the wavenumber a T) can be 
calculated using the procedure outlined at the end of Section 2. The 
functions required in the calculation were evaluated from the series of 
asymptotic expansion of depending on the size of the argument. The 
wavef low equation (2.13) was solved by a finite difference method together 
with an iteration procedure to evaluate the terms involving z derivatives. In 
Figures 4 and 5, we have shown the functions X^(z)^Xj^(z) , 

(i) (l,z) obtained from such a calculation at T = 11,000. We recall that at 
^ yyy ’ 

T = T X =X, = 6,d) (0,z)=u (l,z)=0so that even at about twice the 

col ' yyy yyy 

critical Taylor number the vortices have a significant effect on the waveflow 
problem. We postpone until Section 7 a discussion of the results obtained at 
higher Taylor numbers. The calculation of the finite amplitude Taylor vortex 
beyond T 27,000 was not possible because it is apparently unstable to 
another Taylor vortex mode with wave numbers 2k^. The mode could of course 
be found for T > 28,000 by solving the steady state equations, but such a 
calculation was not carried out. 


6. OBLIQUE WAVES 

The method used to obtain the dispersion relation in Section 2 can be 
extended to deal with Tollmien-Schlichting waves traveling at an angle to the 
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main direction of flow of the Taylor-Gortler vortex. In this case, the 
perturbation velocities and pressure will depend on the slow Tollmien- 
Schlichting cross-stream coordinate ez, as well as the relatively fast 
vortex coordinate z. In the boundary layer, this forces a much larger pres- 
sure gradient in the z-direction which alters the structure of the flow there. 

In the core, then, the perturbation scalings (2.8) remain the same, but 
(2.9) is changed to 


E = h exp(i(aX + 3ez - QT)) h « 1. (6.1) 

This leaves the core flow problem (2.10a-d) unchanged and leads to the same 
matching conditions (2.16 - 2.17). However, because of the z-dependence in 
(6.1), the pressure gradient has a component in the z-direction of O(e^) 

in the boundary layer, as opposed to O(e^) beforehand. This then forces 
the following new scalings in the boundary layer 

y = e^Y, p = P ...]E 

* 2 2 4 5 2 

U = U^[(e XqY,0,0) + (cuq + e Uj^ + . . . , e Vq + e v^ + . . « , e 

3 

e Wj^ ...)E]« 

Substituting the above into the Navier-Stokes equations and linearizing, we 
obtain the following two sets of equations 
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iou« + V 


OY 


+ w 


Oz 


= 0 


and 


i(-« + oXqY)Uq + XqVq + ^oz\q ■ 
l(-fl + aX^Y)wQ = -i8PQ “ Piz 

“ '^0 “ ”0 “ ^ ^ ^ 

Uq. Vq. Wq, >0, Pq = Pq as 

i“Ui + 


“OYY 

OYY 


Y ■*■ “ 


( 6 . 2 ) 


i(-ft + «XqY)Uj + XqV^ + Xq^Yw^ = -iapQ + 
l(-ft + aXQY)Wj = -igp^ - + w^yy 

Uj^ = = Wj^ = 0 Y = 0 

Uj ■>• Aq, Wj^ •>• 0 as Y -»• ®, 

o 2 

Here we have ignored terms of O(h^) and 0(h /e); the latter of the two 
terms makes the extension of this work into the weakly-nonlinear region non- 
trivial. 

Equations (6.2) are the same as (2.10) and (2.20) with (u,v,w) 

replaced by (uq.vqjWq) and ^Pq’^OI ’^0^ replaced by 

z 
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(0, 13Pq + 0). Hence (6.2) can be reduced to the equivalent of (2.25): 

(Pl + iBPQ)^ + iSPo^ “ 

z z 

For periodic solutions we have 




Iz 


■“"oC - ? 


/ 


2ii/k 




dz 


(6.4) 


g 

By adding — x (6.2c) to (6.3b), we find that (6.3) are also the same 

equations as (2.10) and (2.20), but this time (u,v,w) are replaced by 

(^1 Wq.Vj ,wp and (Pq^Pq! replaced by 

z 

(Pod +^) --^Plz’ P2z‘ ^0>* 

a a 


In the case B = 0, equation (6.3) reduce to those for the two dimensional 
disturbance, while (6.2) have the trivial zero solution. For the case 
$ ^ 0 , we obtain the following equation instead of (2.25) 


^P2z ^^Ppz ’*'0^P2z ^^Pp “ ^^^PO ■ ^^Plz 


2 

+ Xq (Io^)Aq 


(6.5) 


(iaXQ)^^\(Ci) 


Substituting for pj^ from (6.4), using the definition of Aq, (2.18), the 
condition for periodic solutions of (6.5) becomes 
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0 


(ct^ I, 


+ -T-) Pr 


<Pi - Pn) Hr 


( 6 . 6 ) 


where Iq, Jq> ^0 Integrals 


^0 ~ 2^ 


2u/k I "^0 

/ e dz 
0 


J 


0 


^ 2ir/k 


2n/k ’*'0 
/ e ^ dz 
0 


(6.7) 


= 


I 

^ /k 


2tt /k 


/ 'J'odz 
0 ^ 


Ai'(5o> 


(iaXQ)'/\(5o) 


dz , 


As in Section 2, we can obtain a similar expression to (6.6) from the boundary 
layer in the upper wall y = 1. Eliminating p^ and p^ between these two 
equations and simplifying leads to the dispersion relation 


1 -tI 


H 


0 


H 


1 


2 I 2 9 9 9 I 

8^/Jq + e^/jj 


( 6 . 8 ) 


where the integrals Ij, Jj, and Hj are defined in a similar manner to (6.7) 
but involving the variables etc., corresponding to the upper 

boundary layer. 

It can be seen that in the case 3 = 0 (6.8) becomes the dispersion 

relation for two dimensional disturbances (2.28). In the case of a two 

dimensional flow with no vortex motion, we obtain Iq = ij = Jq = jj^ = 1 

and Hq = Hj = Ai'(C )/(laX )^^^K(5 ) , so that (6.8) would reduce to 
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(a^ + 6^) = Ai'(0/(iciX)^^^ K(C). 

which is, of course, the usual eigenrelation for three dimensional 
disturbances. We postpone a discussion of the numerical results we have 
obtained for this eigenrelation until the next section. 


7. RESULTS AND DISCUSSION 

We shall concentrate our attention on the effect of longitudinal vortices 
on the growth rate of Tollmien-Schlichtlng waves. Though there is some 
interest in the effect of the vortices on the neutral curve for the Tollmien- 
Schlichting wave, it is the effect of the vortices on the growth rates which 
will be most relevant to the closely related external boundary layer 
problem. In any case our calculations indicate that large amplitude vortices 
have little effect on the neutral configuration whilst even small amplitude 
vortices significantly alter the growth rates in the unstable regime. 

In Figures 6 and 7, we have shown the growth rate and wavenumber of two- 
dimensional Tollmien-Schllchting waves at different values of the Taylor 
number T. The results shown correspond to k = = 3.951 the critical 
wavenumber of linear theory for Taylor-Gortler vortices. The vortices have 
little effect on the neutral frequency and the size of the largest amplifi- 
cation rate. We see in Figure 6 at most frequencies the amplification rate 
increases monotonically with T. For T > 27,000 the Taylor vortex could not 
be calculated because it was apparently unstable to a vortex with twice the 
spanwise wavenumber of the most dangerous mode of linear theory. The 
frequency corresponding to the maximum growth rate increases with T. 
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Moreover, the growth rates beyond the maximum are significantly increased 

for n less than about 100. This result is of particular importance to the 

control of external boundary layers if a similar result holds for such 

flows. Certainly the known similarities between the lower branch structures 

for Polsseuille flow and Blasius flow make that likely, but there are 

difficulties in applying the theory to external flows. The major difficulty 

is surprisingly not the effect of boundary layer growth which can be taken 

• • 

care of as ia Smith (1979b) but the lack of a nonlinear theory for Gortler 
vortices in growing boundary layers. Thus, though our approach of Section 3 
is readily applied to external flows, the absence of any knowledge of even 
weakly nonlinear Gortler vortices at 0(1) wavenumbers prevents us from 

completing such an investigation. 

For external flows, the local Tollmien-Schlichting frequency increases as 

the wave travels downstream and the growth rate adjusts locally. Thus the 

total growth of the disturbance can be found by integrating the growth rate in 

the strearawise direction. In this context the increased growth rates shown in 

Figure 6 to the right of the maximum are possibly significant. As a measure 

of the destabilization produced by the vortices, we can calculate the area 

between the different curves and the “ 0 axis for 20 < ^ < 100. Such a 

calculation shows that for T > 11,000 the area is at least 30% greater than 

that for Poisseulle flow. Thus for external flows which can support 
•• 

Gortler vortices it is possible that their presence might cause the premature 
growth of Tollmien-Schlichting waves. 

We further note that Figure 6 shows that the dependence of on 

becomes increasingly oscillatory when T Increases. We have no physical 
explanation of why that should be the case. The effect of finite amplitude 
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vortices on the Tollmien-Schlicting wavenumber is shown in Figure 7, The 
wavenumber increases monotonically with T, but the rate of increase is very 
small between T = 19000 and T = 27000. 

In Figures 8 and 9, we have shown how the oblique Tollmien-Schlichting 
waves discussed in Section 6 respond to the presence of longitudinal vortices 
at T = 11000. The results are similar to those shown in Figure 4 and again 
suggest that longtitudinal vortices with k « can significantly 
destabilize Tollmien-Schlichting waves. 

We should note that for the channel problem Tollmien-Schlichting 
instabilities might be expected to occur first at finite Reynolds numbers. 
For external flows this is also possible, but there it seems more natural to 
make a high Reynolds number approximation since there would not be a boundary 
layer unless the Reynolds number were large. Thus it might be argued for 
external flows that the most significant linear instability calculation is one 
which calculates the amplification rates between the upper and lower branches 
of the neutral curve. ' Since the motivation for our calculation was to shed 
light on the possible effects of longitudinal vortices on Tollmien-Schlichting 
waves in boundary layers, we feel that a large Reynolds number assumption is 
sensible. We note however that at finite Reynolds the normal and spanwlse 
velocity components of the longitudinal vortex are no longer negligible and 
the z-dependence of the vortex does not become parametric in any region of the 
flow. Thus at finite Reynolds number the computations required would be 
significantly larger than those discussed here. 

Our aim in this work has been to find the effect of finite amplitude 
longitudinal vortex structures on the growth of infinitesimal Tollmien- 
Schlichting waves in curved channel flows. We have Ignored the possibility 
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that the vortices become unstable to time-dependent nonaxisymraetric vortex 
modes of the type which lead to the onset of wavy vortex flows in the Taylor 
problem. We note that Hall ( 1982b) has shown that such disturbances occur in 
external flows over curved walls so this possible mechanism for the onset of a 
time-periodic secondary instability should not be ignored. However, if the 
latter mode does indeed occur in curved channel flows, the question of whether 
it or Tollmien-Schlichting waves are the cause of the secondary instability of 
Taylor-Gortler vortices can only be answered by a nonlinear analysis. 
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CAPTIONS 

Figure 1. The growth rates predicted by the weakly nonlinear theory and the 
corresponding asymptotic results for larger values of ft . The 
full numerical solution for T = 8800 is also shown. 

Figure 2. The functions Ep £ 2 * E^, E^ at different Taylor numbers. 

Figure 3. The functions F^, Fp F 2 , F^, F^ at different Taylor numbers. 

Figure 4. The shear stresses and as functions of z for T = 11,000, 

k = 3.951. 

Figure 5. The functions ()> (-,z), (|i (l,z) for T = 11,000, k = 3.951. 

Figure 6. The wavenumber as a function of ft for several values of 

T^ with 3=0. 

Figure 7. The growth rate as a function of ft for several values of T 

with 3=0. 

Figure 8. The wavenumber as a function of ft for T = 0, 11000 with 

3 = 2. 

Figure 9. The growth rate as a function of ft for T = 0, 11,000 with 


6 = 2 
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